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Abstract 

The paper considers the spectrum of axial perturbations of slowly uni- 
formly rotating general relativistic stars in the framework of Y. Kojima. In a 
first step towards a full analysis only the evolution equations are treated but 
not the constraint. Then it is found that the system is unstable due to a con- 
tinuum of non real eigenvalues. In addition the resolvent of the associated 
generator of time evolution is found to have a special structure which was 
discussed in a previous paper. From this structure it follows the occurrence 
of a continuous part in the spectrum of oscillations at least if the system is 
restricted to a finite space as is done in most numerical investigations. Fi- 
nally, it can be seen that higher order corrections in the rotation frequency 
can qualitatively infiuence the spectrum of the oscillations. As a consequence 
different descriptions of the star which are equivalent to first order could lead 
to different results with respect to the stability of the star. 



1 Introduction 

The discovery |l9j of the instability of r-modes in rotating neutron stars by 
the emission of gravitational waves via the Chandrasekhar-Friedman-Schutz (CFS) 
mechanism |l^ found much interest among astrophysicists. That instability 
might be responsible for slowing down a rapidly rotating, newly-born neutron star 
to rotation rates comparable to the initial period of the Crab pulsar (~19 ms) 
through the emission of current-quadrupole gravitational waves and would explain 
why only slowly-rotating pulsars are associated with supernova remnants 1^, 
Also, while an initially rapidly rotating star spins down, an energy equivalent to 
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roughly 1% of a solar mass would be radiated in the form of gravitational waves, 
making the process an interesting source of detectable gravitational waves Is^] . 



It was soon realized in a large number of studies that many effects work against the 
growth of the r-modc like viscous damping, coupling to a crust, magnetic fields, 
differential rotation and exotic structure in the core of the neutron star. Those 
effects could lead to a significant reduction of the impact of the instability or even 
to its complete suppression. For an account of those studies we refer to the recent 
reviews |||, 

Most of the results have been obtained using a newtonian description of the fluid 
and including radiation reaction effects by the standard multipole formula. Of 
course, such an ad hoc method cannot substitute a fully general relativistic treat- 
ment of the system, but it was believed that it gives at least roughly correct results, 
both, qualitatively and quantitatively. However in a first step towards such a fully 
relativistic treatment it was shown in that the method misses important rel- 
ativistic effects. Working in the low-frequency approximation Kojima could show 
that the frame dragging leads to the occurrence of a continuous part in the spec- 
trum of the oscillations. This is qualitatively different from the newtonian case 
where this spectrum is 'discrete'. Mathematically, Kojima's arguments were not 
conclusive since drawn from an analogy to the equations occurring in the stability 
discussion of non-relativistic rotating ideal fluids |l^ and because his mathe- 
matical reasoning still referred to 'eigenvalues' (which generally don't exist in that 
case) rather than to values from a continuous part of a 'spectrum'. But soon af- 
terwards in | [lO| K. Kokkotas and myself provided a rigorous interpretation along 
with a proof of Kojima's claim. Indications for the continuous spectrum were also 
found in the subsequent numerical investigation . After that still the possibility 
remained that the result was an artefact of the used low-frequency approximation 
which in particular neglects gravitational radiation although numerical results in 
p7t suggested that this is not the case. Also is Kojima's 'master equation' (see 
(]lq)) time independent whereas mathematically it is preferable to start from a time 
dependent equation, because in this way it can be build on the known connection 
between the spectrum of the generator of time evolution and the stability of the 
system For these reasons we consider here Kojima's full equations for 

r-modes which include gravitational radiation but still neglect the coupling to the 
polar modes. Indeed we will meet some surprises, related to the following. 

Due to lack of appropriate exact background solutions of Einstein's field equa- 
tions the background model and its perturbation are expanded simultaneously into 

^But note that depending on the cauation of state still mode solutions can be found in the 
low-frequency approximation [M, B6l M| . 
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powers of the angular velocity of the uniformly rotating star. In particular Ko- 
jima's equations are correct only to first order in Q. Of course, once the evolution 
equations of the perturbations are given, there is no room for neglecting any sec- 
ond nor higher order terms occuring in further computations. The spectrum of 
the oscillations is determined by those equations and depending on that the sys- 
tem is stable or unstable. Now in the calculation it turns out that second order 
corrections in the coefficients of the evolution equations can qualitatively influence 
the spectrum of the oscillations. In particular continuous parts in the spectrum 
(both, stable and unstable) can come and go depending on such corrections. As 
a consequence different descriptions of the star which are equivalent to first order 
could lead to different results with respect to its stability. Hence the decision on 
the stability would have to take into account such corrections which in turn would 
lead to considering a changed operator. In addition judging from the mathemat- 
ical mechanism how this can happen it does not seem likely that this property of 
the equations is going to change in higher orders, which would ultimately question 
the expansion method as a proper means to investigate the stability of a rotating 
relativistic star. With respect to this point further study is necessary, but still the 
results shed some doubts on the appropriateness of the expansion method. 



2 Mathematical Introduction 

Continuous spectra have been found in many cases in the past in the study of 
differentially rotating fluids. ||. The continuous spectrum in these cases was 
again seen for r-modes together with many interesting features such as: the pas- 
sage of low-order r-modes from the discrete part into the continuous part as the 
differential rotation increases; and the presence of low order discrete p-modes in 
the middle of the continuous part in the more rapidly rotating disks |38[ . 

The stars under consideration here have no differential rotation and the existence 
of a continuous part of the spectrum is attributed to the dragging of the inertial 
frames due to general relativity which is an effect not present in an newtonian de- 
scription. 

Mathematically, the study of spectra containing continuous parts requires a higher 
level of mathematical sophistication than usual in astrophysics. Such parts can 
cause instabilities and they cannot be computed by straighforward mode calcula- 
tions. Their occurrence makes it necessary to differentiate between 'eigenvalues' 
(and the corresponding 'modes') and the 'spectrum' of the oscillations. The last 
depends on the introduction of a function space, a topology and the domain of 
definition of a linear operator (namely the generator of time evolution) analogous 
to quantum theory and for this the use of subtle mathematics from functional anal- 
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ysis, operator theory and in particular 'Semigroups of linear operators' is essential. 

Since the system considered here is dissipating (by gravitational radiation) the gen- 
erator of time evolution has complex spectral values and is non self-adjoint. This 
complicates the investigation, because a general spectral theory for such operators 
comparable to that for self-adjoint operators is still far from existing. Indeed here 
there was not (even) found a 'small' self-adjoint part of the generator which would 
have been suitable for applying the usual perturbation methods, n 



3 Kojima's Equations for R-modes 

Since the calculations here are based on the equations of Kojima [p6| which are 
presented in detail there, here we are going only briefly to describe the perturbation 
equations. The star is assumed to be uniformly rotating with angular velocity 
n ~ 0(e) where 



e:=ny/WjM (1) 

is small compared to unity. ^ Here M and R are the mass and the radius of the 
star. Note that we use here and throughout the paper geometrical units c = G = 1. 

The background metric is given by: 

9ab = -e'^dt^ + e^dr^ + r^ {dd^ + sin^ Odcti^) - 2ujr'^ sin^ Odtdcj) , (2) 

where uj ^ 0{e) describes the dragging of the inertial frame. If we include the effects 
of rotation only to order e the fluid is still spherical, because the deformation is of 
order ||2^. The star is described by the standard Tolman-Oppenheimer-Volkov 
(TOV) equations (cf Chapter 23.5 |321) plus an equation for uj 

{jr^vo')' - 167r(p + p)e^jr'^w = , (3) 

where we have defined 

w = Vl-uj (4) 
a prime denotes derivative with respect to r, and 

J = e-(^+'^)/2. (5) 



^ However note that such a method was successfully applied to Kojima's equation ( [l6[ ) for the 
low-frequency approximation [ho||. 

^ The assumption of slow rotation is considered to be a quite robust approximation, because 
the expansion parameter e is usually very small and the fastest spinning known pulsar has e ~ 0.3. 
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In the vacuum outside the star w can be written 



2J 

^ = ^ - 73- ' (6) 

where J is the angular momentum of the star. The function tu, both inside and 
outside the star is a function of r only and continuity of w at the boundary (surface 
of the star, r = R) requires that ro^ = &JR~^. Additionally, w is monotonically 
increasing function of r limited to 

070 < < O, (7) 
where vjq is the value at the center. 

The basic variables for describing r-modes propagating on the background met- 
ric 5^°^ are the functions hoim{t,r),hiim{t,r) {r,6,tp spherical coordinates, t time 
coordinate) defined by expansion into spherical harmonics (imposing the Regge- 
Wheeler gauge) 

he = hot = -^hoimYim,v>/sinO 
hip = /i^pt = 5Z ^oimSinOYi^fi 

l.ra 

he = her = -^hiimYim,,p/sin6 

Lrn 

hr^ = h^r = ^ hiimSinOYirafi , (8) 
l,Tn 

and the fluid perturbation Uim{t,r). Here 

9ab = Qab + hb , (9) 

where hab is the 'small' perturbation. The dependance of the basic variables on 
1,171 will be suppressed in the following. The evolution/constraint equation will be 
written in terms of the following vector- valued variable 

h := I ho 1 (10) 

Then Kojima's equations for pure r-modes {neglecting coupling to polar modes) of 
a slowly and uniformly rotating general relativistic star take the following form 

{I > 2): 

hi = -{Dh)i:=ho- Z (11) 
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ho 



— [2M + kr^ {p - p)] hi + e''-^h[ 



im 
im 



AM 



-A 



ho 



—vje^^r{rZ' - 2h'„) - — [w'e^^r^ + zu {kr^(p + p) - 2re 



A 



(12) 



Z 



iDh)3 

3 [r(A - 2) + 4M + 2kr^{p - p)] hi + -e^'^h'i 



+2— rroe^^Z' + — (Acj - 4e 



-A 



A 

im 



zu) h^ 



+2^ [{2e-^ - kr^{p + p)) m - Alu - e-^ruj'] Z 



+2 



im 

T 



{l + 2e-^)uj' + {2kr^{p + p) - Ar + 4M) 



ZD 



ho 



where A = l{l + 1), k =^ 47r, or in a more compact notation 

'1^-Dh. 
dt 

In addition we have: 

rZ' - 2h'Q+ [2- e^r'^k{p + p)] Z - ire^+^U 

+ — [4M - Ar - 2kr^{p + p)] ho - imrujh'i 



(2 - kr^e^{p + p)) lu + no 



,A + 2 



A 



hi = 



(13) 



(14) 



(15) 



which gives the fluid velocity U in terms of ho, hi and Z. Since U has to vanish 
outside (llq) constrains the data for (04) outside the star. 



3.1 Reminder on Results in the Low- Frequency Approxima- 
tion 

Kojima |2^ investigates the r-modes of the system with low-frequency of the order 
0(e). He finds that the master equation governing those oscillations is given by 



q^+{zu - fi) 



1 







(16) 
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where 



and 



<J> = ^, (17) 



v = ^{l-l)il + 2) , (18) 



l^i- i^^J"^')' = 167r(p + p)e^w , (19) 



2m 



-{a~mn). (20) 



Mainly from its similarity with equations describing plane ideal newtonian fluids 
[ |l^ , [ist he concludes that the spectrum of the oscillations is given by the singular 
values ^ of ( pj| ) inside the star, i.e., by the range 

£(£ + 1) 

wq < fJ- = — — (cT — mQ) < wji . (21) 

In fact in it was proven that it is given by the larger set 

til + I) 

vJo<tJ-^ — ^ (cr - mQ) < n . (22) 
Note that such singular values are not visible in (p^). 

4 The Evolution Equations 

In a first step we deal with the evolution equations only. To formulate a well posed 
initial value problem for the system (|lj) data will be taken from the Hilbert space 

E 

X -.^ Ll{I,j) X Ll{I,j) X Llil,r^j) , /:=(0,oo), (23) 
where j is defined by (^). Further we define an operator A : D{A) X hy 

D{A) := |/I e {C\l,C) X C^{I,C) x C\l,C})r\X : Dh e x} (24) 

*i.e., zeros of the coefficient multiplying the highest order derivatives of the equation 
^For the used notation compare the Conventions in the Appendix. 
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and 



Ah:= Dh , h(£ D{A) . (25) 

Obviously, it can be seen by partial integration that the adjoint operator to A is 
densely-defined ^ and hence that A is closable, i.e., there is a unique 'smallest' 
closed extension of A which is denoted by A in the following. In the following the 
system ( p^ ) is interpreted as the abstract equation 

h{t) ^ -Ah{t) , teR , (26) 

where the dot denotes the ordinary derivative of functions assuming values in X. 
The use of this formulation makes possible the application of the results in the 
field of 'Semigroups of linear operators'. In order that (|26| ) has a unique solution 
for data from the domain of A it has to be proven that A is the generator of 
strongly continuous semigroup (or group). That this is not just a technicality is 
already indicated by the fact that the system ( p^ ) can be seen to have complex 
characteristics if the equation 

r 1 . 

-; ^ wtn < — (27) 

A r~2M f ^ ' 

is violated. Hence in those cases it cannot be expected that A is such a generator. 
For this reason (|2^) is assumed to hold for now on. Note that since 



and because of 



and 



< e~**^Jo ^ iP+P)e t-r (28) 



ujvj ^ —— uja{^ — Wo) for r — !■ (29) 



K r-2M A 



2rn^VlJ 1 , 
u)w ~ tor r ^ oo (30) 



A r-2M A 

that there is a large range of values for the physical parameters where inequality 
( P?! ) is satisfied. But note also that we meet here for the first time a quantity of 
the order O(e^) which restricts the meaningfuUness of (^4|) which itself is correct 
only to the order 0(e). 

^i.e, is defined on a subspace of X which is in addition such that any element of X is a limit 
point of that subspace. 
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If A is the generator of strongly continuous semigroup (or group) -which is to 
expect- its spectrum a{A) ^ tells us about the stability of the solutions of (p6|). 
For instance, if that spectrum contains values from the left half-plane of the com- 



plex plane then there are solutions of ( 26 ) which are exponentially growing. Now 



the whole process of proving that A is such a generator would be greatly simplified 
by the usual perturbation methods if A could be split into a symmetric differential 
operator and into a 'small ' perturbation. This was tried on a diagonalized form of 
the evolution equations (|7^, but unsuccessfully. Indeed the problem there occurred 
in proving the 'smallness' of the the perturbing operator. A further problem in that 
case was that assumptions of usually used theorems giving the asymptotics of solu- 
tions near oo which are needed for the construction of the resolvent of the operator 
turned out to be not satisfied. For these reasons that approach was not pursued any 
further. 

Usually, after the proof that an operator is the generator of a strongly continuous 
semigroup the far more difficult problem occurs in finding physically interesting 
properties of its spectrum and here it might be thought on first sight that this is 
analytically impossible, because of the complicated nature of the evolution equa- 
tions (^4|). But fortunately in |^ there was found a whole class of operators ^ 
which generally have a continuous part in their spectrum and where the occur- 
rence of that part can be concluded from the structure of the resolvent. Indeed in 
the next section the same structure will be found in the resolvent of A. 



5 Construction of the Resolvent of the Generator 

In a first step we try to invert the equations 

{A - a)h = f (31) 

for h G D{A) for any complex number a (= ia in Kojima's notation) and any 
continuous function 

//A 

(32) 

assuming values in C'^ and with compact support on the half line. Any a for which 
there is some / such that ( ^ij ) has not a unique solution is a spectral value. The 

^given by all complex numbers a for which the corresponding map A — a : D{A) X is not 
bijective 

*The operators of that class occurred as a natural generalization of the operators governing 
spheroidal oscillations of adiabatic spherical newtonian stars. 




9 



unique inversion can fail for two reasons. Either there is a non trivial solution of 
the associated homogeneous equation and hence a is an eigenvalue or otherwise 
there is no solution at all for some particular non trivial /. Often the last happens 
not only for some isolated value of a but for a 'continuous set' of values (like a real 
interval, a curve in C or even from an open subset of C) leading to a 'continuous 
part' in the spectrum. For all other values of a the inversion leads to a continuous 
('bounded') linear operator on X. 

Due to the special structure of A (the orders of differentiation inside A vary in 
a special way) these equations can be decoupled leading to a single second differen- 
tial equation for ho. This equation generalizes Kojima's equation (|l^). It is given 

by 



Pi 




(33) 



where 



Pi 

P2 
P3 

Pi 



A 

im r , 



e '^r^ + kr^m {p - 



P)] 



^ ( 



ima 
imuj 



■wr e 



(34) 



91 

<?2 
93 
94 



and 



Auj - 2e~^ruj' ~ vj [ 2kr'^{p + p) 



AM 



A 



m 



^ [2M + kr'^{p - p)] - [uj'e-^r'^ + w {kr^{p + p) - 2re"^)] 



8M 



Auj 



|^(A-2)+ (|-zmw)t7r 



(35) 



91 

92 



-'^r^e-^zuf[ + ^ [cV^e-^ + n7 (fcr^(p + p) - 2re-')] h f2 



' (7 \ r 

-r I - - imu /i + /2 - -/s • 



(36) 
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From ho the functions hi and Z can be computed by 



hi 
Z 

Note that q4 vanishes if only if 



92 
94 



P4_ 
94 



h' 



94 



ho 



/iq + crhi + fi 



muj ± \ m^Lo^ + (A — 2)- 



for some r > and that 
Pi 94 -P3P4 = - 
vanishes if and only if 



1/-A 



1 —ujwre 

A 



17 — -zu 
A 



{n 



1 _ Ijliuwr'^e''' 



(37) 



(38) 



(39) 



(40) 



for some r > 0. Note that the last formula gives up to first order in m exactly the 
values of the continuous spectrum found for (|l^). Also note that the denominator 
in (Eol) is greater than zero, because of 



m 2 

1 —ujmr e 

A 



A r~2M 



(41) 



and condition ( |27|) demanding that the right hand side of the last equation is 
greater than zero^ 

Hence both functions vanish only for purely imaginary a. So the equations are 
non singular for non purely imaginary cr and this case is considered in the follow- 
ing. We denote by Pi the coefhcient of the leading order derivative in (33) and by 
P2 the coefficient of the first order derivative 



Pi 

P2 



P3P4. 




93 
94 



92P4 
94 



(42) 



^Remember that the last condition was imposed to exclude the occurrence of complex charac- 
teristics for the evolution equations dl4h. 
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If hoi and hor (here, ' I ' stands for ' left ' and ' r ' for ' right ') are hnear indepen- 
dent solutions of the homogeneous equation associated with ( p3| ) which are square 
integrable near and near 00, respectively, ho is given by 



ho{r) 



1 

C 



hoi{r) / hor{r')g'i{r')K{r')dr' 

J r 

+ hor{r) f hoi{r')g^{r')K{r')dr' 

JQ 

where C is a non zero constant defined by 



hoi{r)hor{r) - hor{r)hQi{r) = C exp 



P2{r'] 
Piir'Y 



dr' 



and 



K{r) 



Pi{r) 



exp 



^dr' 
Pi{r') 



(43) 



(44) 



(45) 



The lower constant of integration has to be the same in the last two formulas. It 
is kept fixed in the following although its precise value does not enter the formulas 
in any essential way. Note that the inhomogenity in ( p3| ) allows the following 
decomposition into terms containing derivatives of the components of / and terms 
without such derivatives by 



A 



54 



(46) 



where 



g,= '-^[oj 'r^e-^ + vo {kr^p + p) - 2re-^)] h h -.92 • (47) 
A (74 
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Finally from (44),(^) we get by partial integration 



hoir) 



im 1 



hoi(r) / {hoye-^mK)'{r')fi{r')dr 

J r 

hor{r) [ {hoir^e-^vjK)' {r')h{r')dr' 



1 

"c 



1 

"c 



1 

c 



hoiir) / h;,,{r')P^{r')g2{r')K{r')dr' 

Jr 94 

nOO 

Jr 1i 

hoiir) r hor{r')^^^{r')g2{r')dr'+ 

Jr 94 

hor{r) r hoi{r')^^^ir')g2{r')dr' 

Jr 1i 

hoiir) / hor{r')94{r')K{r')dr'+ 

J r 

hor{r) I h,n{r')gAr')K{r')dr' 



Further it follows from (43), and that 



1 

C 



hoiir) / hor{r')g3{r')K{r')dr'+ 

horir) / hoi{r')gAr')K{r')dr' 
Jo 

and hence by (H), (Eq) and partial integration 



^o(^) ^ h^j{r) - -^92 



A Pi 



-/i , 
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where 



iml 



J r 

Kr{r) r (hoir'e-^ujK) ' {r')h{r')dr' 



1 



Ki(r) I KAr')-{r')g2{r')K{r')dr'+ 
94 



Krir) I K,{r'f-^{r')g,{r')K{r')dr' 

94 



+ 



C 



KM horir') 



94 



{r')g2{r')dr'+ 



hM I hoi{r')^P^{r')g2{r')dr' 



94 



Kiir) / W(r-')ff4(r')if(r')d^'+ 
'ior('') / hoi{r')gi{r')K{r')dr' . 



1 



Finally, using ( ^0|) in ( p7|) and some calculation leads to 



im r^e ^vj I imojr 



A 54^1 



and 



where 



Z = H'qj + ahii 



94P1 



2 /i + /2 - ^-^-^ 



— ^ — Ji + J2 - ^•'S 



till ~ n-oi no ■ 

94 94 



(51) 



(52) 



(53) 



(54) 



Note that the structure of these formulas for hi, Z is similar to that of formula (43) 
for ho, with one important difference. Apart from terms containing integrals both 
include an additive term which does not involve integration. They are 



im r^e 1 imojr 



A 94P1 



2 — /i + h ^ -j^fz 



/-A 



94P1 



imujr r 
— ^ — h + /2 - -j^h 



for hi 
for Z . 



(55) 
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Note that both factors multiplying the brackets diverge at the values of a given by 
) . So here we recover again up to first order in m the values of the continuous 
spectrum found for (p^). It is known from |^ that such values become part of 
the spectrum if the operator is considered on a compact interval. |^ For the 
construction of the resolvent there are still needed solutions of the homogeneous 
equation corresponding to (|3^) which are square integrable near r = and others 
which are square integrable near oo. The following Section investigates on their 
existence. 



6 Asymptotics of the Homogeneous Solutions 

The result of the investigation is as follows. In the following we drop the assumption 
that <T is non purely imaginary. There are ^ linearly independent solutions /iqij ^02 
satisfying 

lim r-('+i)/ioiW = 1 
limr-'/if,i(r) = ? + l 

r— >0 

lim r^ho2{r) = 1 

limr'+i/i' (r) = -/ (56) 

and for a different from 

0, ^A(A - 2), ^ [A(A - 2) ± 8mMn] (57) 

there are |^ linearly independent solutions 

ho3{r) - rP^e-'^'-Uiir) 

hoi{r) = r''^eT^''C/2(r) , (58) 

such that 

lim Un{r) = 1 

r — ^00 

lim J7^(r) = . (59) 

r — ^00 

for n — 1,2. Here 71, 72 are solutions of 

7 r^^ —, ^7 - cr^ = (60) 



^''See the proof of Theorem 17 in th 
compactness of integral operators in i 



paper. Basis for this Lemma 2 in the Appendix and the 



and 



^^for e.g., according to the variant of Dunkel's theorem [ ]16[ | given in the Appendix 
^^for e.g., according to the proof of Theorem 1 in paragraph 8 of |24| . 
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such that Re (71) < Re (72), 



CTl7r! 



2jn + O-Q 

for n — 1,2, where 

CTi = {32m^Mn'^a [- (AM + m^nj) a + imMVL{K - 2)] 



(61) 



-2 A \ka - im 



A^ 



a-!^(A-2)17 
A 



CTO 



2cr2 [{Jflm^ + 2MK)a - 2iAmAr 

A[cr-if{A^2)n] 
lem^Ml^^ cr 
A2 



a~if{A~2)n 

Note that the presence of a second order term in (^) diverging near 

a=^(A-2)f^, 



(62) 



(63) 



which is the newtonian frequency for r-modes (to first order in f2) as seen from an 
inertial observer |lOj . Despite of being second order that term becomes arbitrarily 
large near this frequency. 



Further, note that (p7[), (p)6|), (58) imply that the corresponding hn, /112, his, ^14, Z\, 
Z2, Z^, Zi satisfy 

lim r-('+2^/iii(r) = -(cr - imuj{Q)){j(Q)f / {I + 2) 
lim r'-i/ii2(r) = (ff - imuj{Q)){j{0)f / {I - 1) 

— >0 

limr"'Zi(r) = 1 + 1 



lim r'+^Z2(r) ^ 



(64) 



and 



lim r e '^^^h\-s[r) = —^i/a 

r — ^00 

lim r^^^e ''^'7ii4(r) = — 72/c 

r — >cxD 

lim r^~P^e"^^'''Zj.(r) = 2 



lim ri-''^e-''^''Z4(r) 



(65) 
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The outcome for the solutions near is, of course, exactly as expected. As a 
consequence for any a the triple 

hi := *'{hii,hoi,Zi) (66) 

is in X near 0. 



The analysis of the solutions of ( |60| ) uses Theorem (39,1) in |Q which general- 
izes the Routh-Hurwitz theorem to polynomials with complex coefficients. For this 
solutions with Re{j) < 0, Re{'y) > 0, will be referred to as 'stable' and 'unstable', 
respectively. For the analysis define the discrimants Ai and A2 of (pw by 



Ai 



A2 



A 



■{K-2)n 



02 



2A 



(A - 2)f7 



m 
4A2 



(A-2)2f72 



and 



A. 



A2 



A 



■(A -2)17 



(67) 



U2 (^(A - 2)r! - az)] Ai + 



(68) 



Then if, both, Ai 7^ and A2 ^ the number of zeros of ( pO[ ) in the open left 
half-plane of the complex half-plane is given by the sign changes in the sequence 
S' = 1, Ai, A2 and the number of zeros in the open right half-plane is given by the 
number of permanences of sign in this sequence. Note that Ai is > (< 0) outside 
(inside) the circle 



'^2 -^(A- 2)17 



4A2 



(69) 



So we consider the following cases. 



• The first case (corresponding to regions I in Fig.l) considers the values of a 
in the complement of the closed strip M x [0,to(A — 2)ri/A]. Here we have 
Ai > and hence maximal one sign change in the sequence. In particular, 
there is no change in sign near the imaginary axis. So there is maximally one 
stable solution and especially no stable solution near the imaginary axis. 

• The second case (corresponding to region II in Fig 1) considers the values of 
cr in the open strip R x (0, m(A — 2)il/A). Here we have three subcases a), 
b) and c) (corresponding to regions A, B and C, respectively, in Fig 1.) 
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Subcase a) considers those values with Ai > 0. Obviously, this implies 
A2 < 0. As a consequence there is exactly one sign change in the sequence 
and hence there are exactly one stable and one unstable solution. 

Subcase b) considers those values with A2 > 0. This implies Ai < and 
two sign changes in the sequence. Hence in this case there are only stable 
solutions. 

The last subcase c) considers those values satisfying, both, Ai < and 
A2 < 0. Then there is only one sign change and hence exactly one stable and 
one unstable solution. 

Finally, by an application of Rouche's theorem (see e.g. Theorem 3.8 in Chapter 
V of follows that 

• For 

|Re(a)|>^(A-2 + »^) (70) 

there exist, both, a solution with 71 < and with 71 > and hence a stable 
and an unstable solution. 

The asymptotics near 00 of hos, for a near and on the imaginary axis is sur- 
prising. Expected was that in the complement of imaginary axis there are always, 
both, a stable and an unstable solution of ( |60| ) and on the imaginary axis that both 
solutions are purely imaginary. Indeed this would have been the case if the second 
order term in this equation were absent. There one has either exponential growth 
of both solutions (in Regions D) or exponential decay of both solutions (Region 
B). As a consequence each a d B ~ different from the values given in ( p7| ) - is an 
eigenvalue of A. Hence there is a continuum of eigenvalues for A in the open left 
half-plane and hence the evolution given by (|l]), (|l|) and (^ is unstable. 

7 Discussion and Open Problems 

In the previous Section it was found that the spectrum of A contains a continuum of 
unstable eigenvalues leading to an unstable evolution given by (|ll|), (^2|) and (|l3|). 
In addition there was found a continuum of values a (from Regions D in Fig.l - 
which include the spectral values found in the low-frequency approximation as well 
as (^ ) - for which any non trivial solution of the homogeneous equation corre- 
sponding to (|3^) growths exponentially near 00. Hence it is to expect that those 
values are also part of the spectrum of A and hence contribute to the instability. 
Note in particular that their associated growth times displayed in Fig.l can be of 
larger size than those for the found unstable eigenvalues (Region B in Fig.l) which 



18 



Imo 





1.5- 

+ 


D 
+ 


I 




1 ( 




II 


A 


(V 


+/- 


A 




\C\ 


/c 






/y+ 


- 0.2 / 

- + 

' -1- 


\ 0.2 

D 


I 



Re a 



Fig. 1: Sign changes in S for m = 1, M = 1, Q = 1, A = 6. 



suggests that the influence on the evolution dominantes over the influence of those 
eigenvalues. 

A further question is how the constraint (|l5|) - which has not been treated so 
far - is likely to influence the spectrum. The constraint should lead to closed sub- 
space of X being invariant under A. Then the spectrum of the restriction of A to 
that subspace is a subset of the spectrum of A. This should have the effect that 
the eigenvalues in region B (see Fig.l) are becoming discrete instead whereas the 
continuous part of the spectrum in Region D might be left unchanged, because in 
general continuous spectra are less sensitive to such a operation. In any case the 
singular structure of the resolvent at the values given by ( ^o|) will be unchanged 
having the effect that those values remain part of the spectrum at least when the 
system is restricted to finite space as is necessarily done in most numerical investi- 
gations. To decide these questions further study is needed. Note in this connection 
that the occurrence of an unbounded spectrum for A would not be a surprise, 
because for e.g., it is known that a self-adjoint operator has a bounded spectrum 
if and only if it is defined on the whole Hilbert space and in general differential 
operators cannot be defined on the whole of a weighted L^-space. Also would the 
occurrence of an infinitely extended continuous part in the spectrum of oscillations 
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be not surprising since that is quite common for infinitely extended physical sys- 
tems. 



From a physical point of view the main worrying feature of the results is their 
qualitative dependance on second order terms like that in ( pOj) . It is to expect that 
changes in second order to the coefRcents in (|l^) , (|2|) and ( |13| ) influence that term 
which can lead to a qualitative change of the spectrum. As a consequence different 
descriptions of the star which are equivalent to first order could lead to different 
statements on the stability of the star. Hence to decide that stability such correc- 
tions of the coefhcients would have to be taken into account which in turn would 
lead to considering a changed operator. In addition judging from the mathematical 
mechanism how the coefficients of the evolution equations influence the spectrum 
- namely through the asymptotics of the homogeneous solutions of near oo - 
it might be suspected that this property of the equations does not change in higher 
orders, which would ultimately question the expansion method as a proper means 
to investigate the stability of a rotating relativistic star. Also with respect to these 
points further study is necessary, but still the results raise first doubts whether the 
slow rotation approximation is appropriate for this purpose. 



A final interesting question to ask is whether a numerical investigation of the evolu- 
tion given by (^2|) and (13) is capable of detecting the computed instabilities. 
This seems unlikely, because in that process space has to be 'cutoff' near the singu- 
lar points r = 0, cx) and suitable local boundary conditions have to be posed at that 
the endpoints. But such a system is qualitatively different from the infinite system 
since for instance there the asymptotics of the homogeneous solutions (^Tj) near oo 
does not play a role. Note that - independent from the used boundary conditions 
- for such a system the continuum of values given by the restriction of (|4^) to the 
chosen interval are part of the spectrum |^ of the corresponding operator as also is 
also indicated by the numerical investigation p7[|. 
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8 Appendix 
8.1 Conventions 

The symbols N, M, C denote the natural numbers (including zero), all real numbers 
and all complex numbers, respectively. 

To ease understanding we follow common abuse of notation and don't differen- 
tiate between coordinate maps and coordinates. For instance, interchangeably r 
will denote some real number greater than or the coordinate projection onto 
the open interval / := (0, oo). The definition used will be clear from the context. 
In addition we assume composition of maps (which includes addition, multipli- 
cation etc.) always to be maximally defined. So for instance the addition of two 
maps (if at all possible) is defined on the intersection of the corresponding domains. 

For each /c G N \ {0}, n G N \ {0} and each non-trivial open subset M of M" 
the symbol (M, C) denotes the linear space of /c-times continuously differen- 
tiable complex- valued functions on M. Further Cq{M, C) denotes the subspace of 
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C'^ (M, C) consisting of those functions which in addition have a compact support 
in f2. 



Throughout the paper Lebesgue integration theory is used in the formulation of 
p5| | . Compare also Chapter III in and Appendix A in |Q . To improve read- 
ability we follow common usage and don't differentiate between an almost every- 
where (with respect to the chosen measure) defined function / and the associated 
equivalence class (consisting of all almost everywhere defined functions which differ 
from / only on a set of measure zero). In this sense Lq {M,p), where p is some 
strictly positive real- valued continuous function on M, denotes the Hilbert space 
of complex- valued, square integrable (with respect to the measure pd"'x) functions 
on M. The scalar product < | > on {M, p) is defined by 

</l5>:= / rgpdTx , (71) 
J hi 

for all f,gE {M, p), where * denotes complex conjugation on C. It is a standard 
result of functional analysis that Cq (M, C) is dense in ( Af , p) . 

Finally, throughout the paper standard results and nomenclature of operator the- 
ory is used. For this compare standard textbooks on Functional analysis, e.g.,|Q 
VoL I, 

8.2 Auxiliary Theorems 

The variant of the theorem of Dunkel (compare also |[ ^ ) used in Section 
4 is the following. 

Theorem 1 ; Let n e A^; a, to e i? with a < to; p, G N; := 1 for p = and 
o^fi '■— l^ for p ^ 0. In addition let Ao be a diagonalizable complex n x n matrix 
and e'l, . . . , e'^ be a basis of C" consisting of eigenvectors of Ao . Further, for each 
j G {\, ■ ■ ■ ,n} let \j be the eigenvalue corresponding to e^- and Pj be the matrix 
representing the projection of C" onto C.e'^ with respect to the canonical basis of 
C". Finally, let Ai be a continuous map from {a, to) into the complex n x n ma- 
trices M{n X n, C) for which there is a number c G (a, to) such that the restriction 
of Aijk to [c, to) *s Lebesgue integrable for each j, fc £ 1, ...,n. 

Then there is a map R : (a, to) M{n x n,C) with limt^oRjk{t) — for 
each j,k G 1, . . . , n and such that u : (a, to) — > M{n x n, C) defined by 

r EUito t)--^ ■ {E + R{t)) ■ P, forp = 
^ ^ I EU e^P(^.(*o - ■ {E + R{t)) ■ P, forp^O ^'^> 
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for all t £ (a, to) (where E is the n x n unit matrix), maps into the invertible nx n 
matrices and satisfies 

for each t G {a, to). 



Lemma 2 Let X be a non-trivial Hilbert space with scalar product < | > and 
induced norm \\ \\ and let be A : D{A) ^ X a densely- defined, linear and closable 
operator in X. Further let be a G C such that A — a is injective. Finally, let be 
fi £ <C* and let be rio,rii, ... a sequence of elements of Ran{A — a) for which there 
is an e £ (0, oo) such that \\ri^\\ > s for every v £ N and 



lim 



{A~a) =0. (74) 



Then ct + /i ^ is in the spectrum of the closure A of A. 

Proof: The proof is indirect. Assume otherwise that ct + /Lt~^ is not part of the 
spectrum of A. Then A — (fT + /i~^) is in particular bijective with a bounded inverse 
(A - (cr + ^-i))" . Then we have for every G N 



{A-{<j + fi-^)) {A - a)-' v.=V.- M"' (A - a)-' rj, (75) 
and hence also 

{A ^ {a + {A - (j)-'^ 7]^ = //(i- (CT + /i-i))"%^ - 

n{A-ay\,. (76) 

Using this it follows from (Q) and the continuity of (A — {a + /i^^)) ^ that 

Jirn |(i- ((T + ^-i))"^ {A - ay'^ T]^ - H {A - {cr + ^i-^)y\^} = 
-II lim (A - cr)"S<^ = (77) 

V — >oo 

and hence by ( [f^ that 

lim f]^ = . (78) 

u — *oo 

The last is in contradiction to the assumption there is an e S (0, oo) such that 
WVvW ^ ^ for G N. Hence the Lemma is proven. □ 



25 



